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I | We investigate parity-violating observables in the np system, including the 

longitudinal asymmetry and neutron-spin rotation in np elastic scattering, 
^ ■ the photon asymmetry in ftp radiative capture, and the asymmetries in 

■ deuteron photo-disintegration d(^y, n)p in the threshold region and electro- 
disintegration d(e, e')np in quasi-elastic kinematics. To have an estimate 
of the model dependence for the various predictions, a number of different, 
latest-generation strong-interaction potentials — Argonne vi$, Bonn 2000, and 

i -^h ' Nijmegen I — are used in combination with a weak-interaction potential con- 

sisting of 7T-, and (j-meson exchanges — the model known as DDH. The 
complete bound and scattering problems in the presence of parity-conserving, 
including electromagnetic, and parity-violating potentials is solved in both 
configuration and momentum space. The issue of electromagnetic current 
' conservation is examined carefully. We find large cancellations between the 

■ asymmetries induced by the parity-violating interactions and those arising 
from the associated pion-exchange currents. In the ftp capture, the model 
dependence is nevertheless quite small, because of constraints arising through 
the Siegert evaluation of the relevant E\ matrix elements. In quasi-elastic 
electron scattering these processes are found to be insignificant compared to 
the asymmetry produced by ^-Z interference on individual nucleons. These 
two experiments, then, provide clean probes of different aspects of weak- 
interaction physics associated with parity violation in the np system. Finally, 
we find that the neutron spin rotation in ftp elastic scattering and asymmetry 
in deuteron disintegration by circularly-polarized photons exhibit significant 
sensitivity both to the values used for the weak vector-meson couplings in 
the DDH model and to the input strong-interaction potential adopted in the 
calculation. This reinforces the conclusion that these short-ranged meson cou- 
plings are not in themselves physical observables, rather the parity- violating 
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mixings are the physically relevant parameters. 
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I. INTRODUCTION 



A new generation of experiments have recently been completed, or are presently under 
way or in their planning phase to study the effects of parity- violating (PV) interactions in pp 
elastic scattering [1], np radiative capture [2] and deuteron electro-disintegration [3] at low 
energies. There is also considerable interest in determining the extent to which hadronic weak 
interactions can affect the longitudinal asymmetry measured by the SAMPLE collaboration 
in quasi-elastic scattering of polarized electrons off the deuteron [4], and therefore influence 
the extraction from these data (and those on the proton [5]) of the nucleon's strange magnetic 
and axial- vector form factors at four-momentum transfers squared of 0.04 and 0.09 (GeV/c) 2 . 

The present is the third in a series of papers dealing with the theoretical investigation of 
PV interaction effects in two-nucleon systems. The first [6] was devoted to pp elastic scatter- 
ing, and presented a calculation of the longitudinal asymmetry induced by PV interactions 
in the lab-energy range 0-350 MeV. The second [7] provided a rather cursory account of a 
study of the PV asymmetries in np radiative capture at thermal neutron energies and in 
deuteron electro-disintegration at quasi-elastic kinematics. This work further extends that 
of Ref. [7] by investigating the neutron spin rotation at zero energy and the longitudinal 
asymmetry in np elastic scattering at lab energies between and 350 MeV, and the photon- 
helicity dependence of the d(j, n)p cross-section from threshold up to energies of 20 MeV. 
It also provides a thorough analysis of the results already presented in Ref. [7]. 

We adopt the PV potential developed by Desplanques et al. [8] over twenty years ago, the 
so-called DDH model. In the np sector, it is conveniently parameterized in terms of 7T-, p-, 
and cu-meson exchanges. In Ref. [8] the pion and vector-meson weak coupling constants were 
estimated within a quark model approach incorporating symmetry techniques like SU(6)vk 
and current algebra requirements. Due to limitations inherent to such an analysis, however, 
the coupling constants so determined had rather wide ranges of allowed values. 

Our prime motivations are to develop a systematic and consistent framework for study- 
ing PV observables in the few-nucleon systems, where accurate microscopic calculations are 
feasible, and to use available and forthcoming experimental data on these observables to con- 
strain the strengths of the short- and long-range parts of the two-nucleon weak interaction. 
Indeed, in Ref. [6] we showed, for the case of the longitudinal asymmetry measured in pp 
elastic scattering, how available experimental data provide strong constraints on allowable 
combinations of p- and cj-meson weak coupling constants. 

The remainder of the present paper is organized as follows. In Sec. II the PV potential 
as well as the parity-conserving strong-interaction potentials used in this work are briefly 
discussed, while in Sec. Ill the model for the nuclear electro-weak currents is described, 
including the electromagnetic two-body terms induced by the presence of PV interactions. 
In Sec. IV a self-consistent treatment of the np bound- and scattering-state problems in 
both configuration and momentum spaces is provided, patterned after that of Ref. [6], and 
in Sec. V explicit expressions are derived for the longitudinal asymmetry and spin rotation 
in np elastic scattering, the photon asymmetry in np radiative capture, and the asymmetries 
in deuteron photo- and electro-disintegration. In Sec. VI the techniques used to calculate 
the PV observables are briefly reviewed, while in Sec. VII a fairly detailed analysis of the 
results is offered. Finally, Sec. VIII contains some concluding remarks. 
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II. PARITY-CONSERVING AND PARITY- VIOLATING POTENTIALS 



The parity-conserving (PC), strong-interaction potentials used in the present work are 
the Argonne u 18 (AV18) [9], Nijmegen I (NIJM-I) [10], and CD-Bonn (BONN) [11] mod- 
els. They were discussed in Ref. [6] in connection with the calculation of the longitudinal 
asymmetry in pp elastic scattering. Here, we briefly summarize a few salient points. 

The AV18 and NIJM-I potentials were fitted to the Nijmegen database of 1992 [12,13], 
consisting of 1787 pp and 2514 np scattering data, and both produced y 2 per datum close 
to one. The latest version of the charge-dependent Bonn potential, however, has been fit to 
the 1999 database, consisting of 2932 pp and 3058 np data, for which it gives x 2 P er datum 
of 1.01 and 1.02, respectively [11]. The substantial increase in the number of pp data is due 
to the development of novel experimental techniques — internally polarized gas targets and 
stored, cooled beams. Indeed, using this technology, IUCF has produced a large number 
of pp spin-correlation parameters of very high precision, see for example Ref. [14]. It is 
worth noting that the AV18 potential, as an example, fits the post- 1992 and both pre- and 
post-1992 pp (np) data with x 2 's of 1.74 (1.02) and 1.35 (1.07), respectively [11]. Therefore, 
while the quality of their fits (to the pp data) has deteriorated somewhat in regard to the 
extended 1999-database, the AV18 and NIJM-I models can still be considered "realistic". 

As already mentioned in Sec. I, the form of the parity-violating (PV) weak-interaction 
potential was derived in Ref. [8] — the DDH model. In the isospin space of the np pair it is 
expressed as 

v ?y T = (T', M' T = | v py | T, Mt — 0) , (2.1) 
where T, T'—O or 1. The diagonal and off-diagonal terms are then obtained as 

m Q (l + K a )Y'(m a r)(rri x cr 2 ) • r 

(2.2) 
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(2.3) 



and VqX=v\q . In the equations above the relative position and momentum are defined as 
r = ri — r 2 and p = (pi — p 2 )/2, respectively, [...,...], denotes the anticommutator, and m, 
m n , m p , and m w are the proton, pion, p- meson, w-meson masses, respectively. The Yukawa 
function Y(x a ), suitably modified by the inclusion of monopole form factors, is given by 



Y(x a ) = 



^-(A a /m a )x 



Xr 



1 + 



lAo 

2 m r 



(2.4) 



where Note that Y'(x) denotes dY(x)/dx, and that the terms proportional to 

Y'(x) in Eqs. (2.2) and (2.3) are usually written in the form of a commutator, since 
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i [p , Y(m a r)}_ = m a Y'(m a r) r . (2.5) 

Finally, the values for the strong-interaction n- meson pseudoscalar coupling constant g n , 
and p- and u;-meson vector and tensor coupling constants g a and n a , as well as for the 
cutoff parameters A a , are taken from the BONN model [11], and are listed in Table I. The 
weak-interaction vector-meson coupling constants h n p v and h™ p correspond to the following 
combinations of DDH parameters 

h? = (4T-3)h° p ~h 2 p , (2.6) 
K P = ^ , (2.7) 

where T—0,1. The values for these and for h n , h p , and are listed in Table I. Note that we 
have taken the linear combination of p- and u;-meson weak coupling constants corresponding 
to pp elastic scattering from the earlier analysis [6] of these experiments. The remaining 
couplings are the "best value" estimates, suggested in Ref . [8] . 

In order to study the sensitivity of the calculated PV observables to variations in the 
weak coupling constants, we also consider in the present work a v PV corresponding to the 
original "best value" estimates for these from Ref. [8], see Table II. 



III. ELECTROMAGNETIC AND NEUTRAL WEAK CURRENT OPERATORS 

The electromagnetic current and charge operators, respectively j and p, are expanded 
into a sum of one- and two-body terms 

j(q) = £j*(q) + £ju(q) , (3.i) 

i i<j 

and similarly for p(q). The one-body terms have the standard expressions obtained from a 
non-relativistic reduction of the covariant single- nucleon current [15]. The two-body charge 
operators are those derived in Ref. [16]; they only enter in the calculation of the asymmetry 
in the deuteron electro-disintegration at quasi-elastic kinematics, and will not be discussed 
further here. 

The two-body currents have terms associated with the parity-conserving (PC) and parity- 
violating (PV) components of the interaction, respectively j PC and j PV . The operators j PC 
were derived explicitly in Ref. [17], and a complete listing of those relative to the Argonne Vi$ 
(AV18) interaction [9] has been given most recently in Ref. [15]. Only the two-body currents 
associated with tc- and p-exchange are retained in the case of the Bonn (BONN) [11] and 
Nijmegen-I (NIJM-I) [10] interactions. 

In addition to these, the purely transverse two-body currents associated with the ex- 
citation of A-isobars and the prc-f and urc-f mechanisms are included in all calculations. 
Again explicit expressions for these operators can be found in Ref. [15]. Note, however, that 
the A-isobar degrees of freedom are treated in perturbation theory rather than with the 
transition-correlation-operator method [18], and that only effects due to single A excitation 
are considered, according to Eqs. (2.15) and (3.4) in Ref. [15]. 

Before moving on to a discussion of the PV currents, we briefly review, for later reference, 
the question of conservation of the electromagnetic current for the case of the AV18. As 
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pointed out in Ref. [15], the currents from its vq part (specifically, its isospin-dependent 
central, spin-spin, and tensor components) are strictly conserved. In a one-boson-exchange 
model, which the AV18 is not, these interaction components arise from tt- and p-exchange. 

The currents from the AV18 momentum-dependent (p-dependent) components — the 
spin-orbit, L 2 , and quadratic-spin-orbit terms — are also included. In Ref. [19] and later pa- 
pers, the currents from the spin-orbit term were derived by generalizing the procedure used 
to obtain the v$ currents. It was assumed that the isospin-independent (isospin-dependent) 
central and spin-orbit interactions were due to a and uj exchanges (p exchange), and the 
associated two-body currents were constructed by considering corresponding AA-pair dia- 
grams involving these meson exchanges. The currents from the L 2 and quadratic-spin-orbit 
interactions were obtained, instead, by minimal substitution [15,17]. 

The currents from the p-dependent interactions are strictly not conserved, as one can 
easily surmise by considering their commutator with the charge density operator. For ex- 
ample, in the case of the isospin-dependent L 2 and (L ■ S) 2 interactions, this commutator 
requires the presence of currents with the isospin structure (t^ x Tj) z , which cannot be 
generated by minimal substitution [17]. 

We will return to this issue in Sec. VII B. Here, we only want to emphasize that the 
currents from the p-dependent terms in the AV18 are short-ranged. Their contributions to 
isovector observables, such as, for example, the magnetic form factors of the trinucleons [20], 
are found to be numerically much smaller than those due to the leading t> 6 currents. These 
currents also lead to small, although non- negligible, corrections to isoscalar observables, such 
as the deuteron magnetic moment and B(q) structure function [21]. However, in the case of 
the PV asymmetry in the np radiative capture at thermal energies under consideration in 
the present study, they will turn out to play an important role (see Sec. VII B). 



A. Parity-violating currents 

The DDH PV interaction [8] is parameterized in terms of 7T-, p-, and c^-meson exchanges. 
The meson-nucleon phenomenological Lagrangian densities have been given most recently 
in Ref. [22]. We adopt here the notation and conventions of that work, except that we use 
pseudo-vector coupling for the ttNN interaction Lagrangian, i.e. 

LI C nn = -—Ni,1»tN ■ , (3.2) 

with f 7T /m n =g 7T /(2m). The resulting 'ynNN coupling is given by 

L™ NN = -e^N^Mr x *)zNA" , (3.3) 

and the •ynNN current is then obtained from the Feynman amplitude in Fig. 2, panel (a). 
The complete PV 7r-exchange current is derived from a non-relativistic reduction of both 
amplitudes in Fig. 2, and to leading order reads: 



• PV Or lr \ _ _ f n ^ n I _ \ 

Jn,ij\ K ii K j) — 7k V T * ' T 3 ~ r 2,i r 2J/ 



h 2 — h 2 



[v n (kj) - v n (ki)] (<Ti ■ ki - (t j ■ kj) 



(3.4) 
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where v n (k) is defined as 



v w (k) = 



'A 2 



mz 



k 2 + Al k 2 + ml ' 



(3.5) 



kj=p- — pi is the fractional momentum delivered to nucleon % (with these definitions q=kj + 
kj), and is a short-range cutoff. In the limit of point-like couplings, the current above 
is identical to that listed in Eqs. ( A7a) and ( A9a) of Ref. [22] . iViV-pair terms arising from 
the photon coupling to the nucleon line containing a PV nNN vertex do not contribute to 
leading order. Lastly, we note that jjj! v satisfies current conservation with the PV 7r-exchange 
interaction, given in momentum space by 



.pv, 



kj, kf) — i 



since 



^(ki-q,^), 



Pi 



+ 



V2m, 
Hk l ,k,-q),P 

fir h n 



Ti X Tj) V w (kj) (Ti kj + i # j , 



(3.6) 



V2 



(Ti ■ Tj - T z ,iT zJ ) [V n (ki) ((Ti + (Tj) -ki+i^ j] , (3.7) 



which is easily seen to be the same as (kj + kj) • j^ v (kj,kj). Here Pi denotes the isospin 
projection operator 



l + r ; 



2,1 



(3.8) 



In the present work, the PV currents induced by p- and cj-meson exchanges have been 
neglected, since, due their short-range character, the associated contributions are expected 
to be tiny. We note, however, that in Eq. (A5) of Ref. [22] a ^pNN contact term, originating 
from gauging the pNN tensor coupling, has been ignored. It is given by 



^ p nn = e 9 -^Na, v ( T x pP) z NA v , 



(3.9) 



and leads, in leading order, to an additional term in Eq. (A7c) of Ref. [22] of the form 



2 m 



(Ti X Tj) z V p (kj) 



(T i x a , + 



J n (Tj x k,- 



(3.10) 



with v p (k) defined similarly as in Eq. (3.5). The term above, when combined with that having 
the same structure in the second line of Eq. (A7c), generates a contribution proportional to 
1 + k p which, in view of the large value of the pNN tensor coupling constant (k p —6.6), is 
expected to be dominant in the PV p-exchange current. 

Finally, there is a PV one-body current originating from the nucleon's anapole moment. 
It can be derived, for example, by considering pion-loop diagrams where one of the vertices 
involves a PV nNN coupling; it has the structure, to leading order, given by 



jf V - 



q) 



Q 2 



2 m 2 



a S (ql)+a v (ql)r z ,i 



(Tie 



iq-ri 



(3.11) 
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where q 2 is the four-momentum transfer, q 2 —cu 2 — q 2 , and the isoscalar and isovector anapole 
form factors are normalized as 

a s ' v (0) = -3^ a sy (312) 

with a s —1.6 and a v =0A from a calculation of pion-loop contributions [23]. More recent 
estimates of the nucleon anapole form factors predict [24-26] somewhat different values for 
a s ' v . A complete treatment would require estimates of short-distance contributions [27] 
and electroweak radiative corrections. Thus, in view of the uncertainties in the quantitative 
estimate of these effects, we will continue to use the values above in the present study (see 
also Sec. VII D). Note that jf v vanishes for real-photon transitions. 



B. Neutral weak currents 

In the Standard Model the vector part of the neutral weak current, j°' a , is related to 
the isoscalar (S) and isovector (V) components of the electromagnetic current, denoted 
respectively as jg' a and j v ' <T , via 

j°' a = -2sm 2 6 w3 r + (1 - 2sin 2 M ft" , (3.13) 

where 9 W is the Weinberg angle, and therefore the associated one- and two-body weak charge 
and current operators are easily obtained from those given in Sec III. The axial charge and 
current operators too have one- and two-body terms. Only the axial current 

j 5 (q) = £j?(q) + £j?,-(q) (3-i4) 

i i<j 

is needed in the present work. The one- and two-body operators are essentially those listed 
in Ref. [28], except for obvious changes in the isospin structure having to do with the fact 
that we are dealing here with neutral rather than charge-raising/lowering weak currents, 
and for the inclusion of nucleon and NA axial form factors — the parameterization adopted 
for these is given in Ref. [29]. Note that in Ref. [31] the relativistic corrections in jf and two- 
body axial currents were neglected, in line with the expectation, confirmed in the present 
study, that the associated contributions were small. 

Finally, the neutral weak currents given above are at tree-level, electro-weak radiative 
corrections as well as strange-quark contributions to the vector and axial- vector currents [30] 
have been ignored. These effects have been taken into account in recent calculations of the 
longitudinal asymmetry in d(e,e')pn at quasi-elastic kinematics [32], however, they will not 
be discussed further here. 

IV. FORMALISM 

In this section we discuss the np scattering- and bound-state problems in the presence 
of a potential v given by 

v = v FC + v FV ? 
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where v PC and v PV denote the parity-conserving (PC) and parity-violating (PV) components 
induced by the strong (including electromagnetic) and weak interactions, respectively. The 
formalism and notation are similar to those developed in Ref. [6]. 



A. Partial-wave expansions of scattering state, T- and 5-matrices 

The Lippmann-Schwinger equation for the NN scattering state | p, SMs^T)^, where 
p is the relative momentum, S, M s , T, and M T =0 specify the pair spin, spin-projection, 
isospin, and isospin-projection states (note that the label M T =0 is unnecessary, since v is 
diagonal in M T ), can be written as [33] 

\p,SM s ,T)W=\p,SM s ,T)o + - j. ,. v \p,SM s ,T)W , (4.2) 

hi — tig ± ie 

where H is the free Hamiltonian, and | . . .) are the eigenstates of H , namely plane waves, 
P ,5M Sl T(r) = (r | p,SM s ,T) 



_L [ e iP-r _ (_)5 + r e -ip.r 



S T 

XmsVmt^ 



= 4ttV2 £ i L tLST J L ( P r) [Z J L ^ s (p)Y y%j(*)rfi ■ (4.3) 

JMjL 

Here ji(pr) denotes the regular spherical Bessel function, and the following definitions have 
been introduced: 

ZtsM s (i>) = Y,( LM l, SM S | JMj) Y LML (p) , (4.4) 

M L 



e LS T = - 



The factor €lst ensures that the plane waves are properly antisymmetrized. 
The T-matrix corresponding to the potential v is defined as [33] 

T(p',S'M' s ,T';P,SM s ,T) = (p', S'M' S ,T' \ v | p,SM s ,T)^ . (4.6) 

Insertion of the plane-wave states | p, SM$, T% into the right-hand-side of the Lippmann- 
Schwinger equation leads to 

|p, SM s ,T)^ = \p, SMs,T) 

ST f dp ' 1 In' q>m> T>\ T(p',S'M' s ,T';P,SM s ,T) 

from which the partial wave expansion of the scattering state is easily obtained by first 
noting that the potential, and hence the T-matrix, can be expanded as 



o(p', S'M' S ,T' | v | p,SMs,T) = 2(4tt) 2 ]T Y,*vs>t> e LST Z J L ^ M ,{p') 

JMj LL> 

[Z J lsMp)V vJ l'S'T',lst(p';p) , (4-8) 
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with 

i>s>T>,LST(p';p) = ^'J drj L ip'r)y^ s } jV pv(r)^y^ j3L (pr) . (4.9) 



v, 



After insertion of the corresponding expansion for the T- matrix into Eq. (4.7) and a number 
of standard manipulations, the scattering-state wave function can be written as 

VfiU, r (r) = W2 ]T £ ^ L '^s^e LST [Z J L ^ s m* 

JMj LL'S'T' 

^^ feWj , (4.10) 

with 



= ft^M + |jf W) 8 ,^ ) + fc a;>) , (4.1D 

where the label a (a') stands for the set of quantum numbers LST (L'S'T'). The (complex) 
radial wave function w(r) behaves in the asymptotic region r — > oo as 



it; 



6 al , a h { ;)(pr) + h { l)(pr)S J af>a (p)\ , (4.12) 



r 2 

where the on-shell (p' = p) ^-matrix has been introduced, 

Sa'M = - 4i/ipT Q J , Q (p;p) , (4.13) 

and the functions h^ 1,2 \pr) are defined in terms of the regular and irregular (ul) spherical 
Bessel functions as 

h { L' 2 \pr) = j L (pr) ±in L (pr) . (4.14) 



B. Schrodinger equation, phase-shifts, mixing angles, and the scattering amplitude 

The coupled-channel Schrodinger equations for the radial wave functions w(r) read: 

wZ, a (r;p) + J2 rv *'A r ) l w L( r '>p) = > ( 4 - 15 ) 

3 1 

with 

vijr) = i L ~ L ' 2/i / dn V p v(r) rf, y% , (4.16) 

where, because of time reversal invariance, the matrix v 3 a , a can be shown to be real and 
symmetric (this is the reason for the somewhat unconventional phase factor in Eq. (4.16); in 
order to maintain symmetry for both the v PC - and v PV -matrices, and hence the S'-matrix, the 
states used here differ by a factor i L from those usually used in nucleon-nucleon scattering 



dr s 



+ L^l)_ p2 
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analyses). The asymptotic behavior of the u>(r)'s is given in Eq. (4.12), while explicit 
expressions for the radial functions v^F^(r) can be found in Ref. [6] — those associated with 
v PC are well known. 

There are two coupled channels for J=0, and four coupled channels for J > 1. The 
situation is summarized in Table III. Again because of the invariance under time-inversion 
transformations of t> PC + v PV , the S'-matrix is symmetric (apart from also being unitary), 
and can therefore be written as [33] 

S J = U T S&U, (4.17) 

where U is a real orthogonal matrix, and Sp is a diagonal matrix of the form 

Si, a >, a = 6 a ,, a e 2iS « . (4.18) 

Here 5^ is the (real) phase-shift in channel a, which is function of the energy E with 
p = y/2fxE. The mixing matrix U can be written as 

U = U {12) J = , (4.19) 

= n u m j ^ 1 . ( 4 - 2 °) 

i<i<j<4 

where is the 2 x 2 or 4 x 4 orthogonal matrix, that includes the coupling between 
channels % and j only, for example 





COS 





sin e{ 3 


0" 




" 





1 


0" 


1/(13) = 




-sin e{ 3 


1 






cos e{ 3 






-1+^3 




-1 

























1 . 

















Note that no coupling is allowed between channels 3 and 4 in the notation of Table III, and 
hence = 1. Thus, for J=0 there are two phase-shifts and a mixing angle, while for 

J > 1 there are four phase-shifts and five mixing angles. Of course, since \v PY \ |f PC |, the 
mixing angles induced by t> PV are <C 1, a fact already exploited in the last expression 
above for U. Given the channel ordering in Table III, Table IV specifies which of the channel 
mixings are induced by v PC and which by w PV . 

The reality of the potential matrix elements v^, a {r) makes it possible to construct real 
solutions of the Schrodinger equation (4.15). The problem is reduced to determining the 
relation between these solutions and the complex u>(r)'s functions. Using Eq. (4.17) and 
U T U = 1, the iy(r)'s can be expressed in the asymptotic region as 

r p 1 

Tn . xJ sin[pr - L' n/2 + 8i\ 
= E ( u ^e 1 ^ — PA U Pa . 4.21 

The expression above is real apart from the exp(i<5^). To eliminate this factor, the following 
linear combinations of the u>(r)'s are introduced 
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p r 



(U rr ) a 'a [cos 5^ j L '(pr) - sin <J^n L '(pr) 



(4.22) 



and the u(r)'s are then the sought real solutions of Eq. (4.15). 

The asymptotic behavior of the u(r)'s can now be read off from Eq. (4.22) once the 
[/-matrices above have been constructed. The latter can be written, up to linear terms in 



the "small" mixing angles induced by t> PV , as 



U 



-12 



-12 



J = 



u = 



cose 



12 



-sine 



12 



sine 

cose 



— e 



-14 



— e 



12 

,/ 
23 



-24 



e l3 COSe l2 + e 23 Sme l2 

- e l3 Sme l2 + e 23 COSe 12 
1 





e{ 4 cose{ 2 + e^sine^ 

o 
l 



j > l 



Inverting the first line of Eq. (4.22), 



w:. 



Uf3 a , 



(4.23) 



and inserting the resulting expressions into Eq. (4.15) leads to the Schrodinger equations 
satisfied by the (real) functions u(r). They are identical to those of Eq. (4.15), but for the 
u>(r)'s being replaced by the tt(r)'s. These equations are then solved by standard numerical 
techniques. Note that: i) v J oia = d^ c , since the diagonal matrix elements of t> PV vanish 



because of parity selection rules; ii) terms of the type r X r )up a (r) / r involving the 
product of a PV potential matrix element with a v PV -induced wave function are neglected. 

Finally, the physical amplitude for np elastic scattering from an initial state with spin 
projections m n , m p to a final state with spin projections m' n , m' p is given by 



(m>; I M | m n m p ) = - £ (-) T+T (-m' n , -m' p \ S'M' s )(-m n , -m p \ SM S ) 



n p I o I \ rj cy 

" s'm' s t',sm s t " 
Ms'm' s t',sm s t(E,9) , 

where the amplitude M is related to the T-matrix defined in Eq. (4.6) via 
M s , M , T , tSMsT (E,6) = -^-T(p',S'M' s ,T';pZ,SM s ,T) , 



(4.24) 



(4.25) 



and the factor (— ) T+T '/2 comes from the Clebsch-Gordan coefficients combining the neutron 
and proton states to total initial (final) isospin T (T') . Note that the direction of the initial 
momentum p has been taken to define the spin quantization axis (the z-axis), 9 is the 
angle between p and p', the direction of the final momentum, and the energy E = p 2 / (2/i) 
(— P' 2 /(2/u))- Using the expansion of the T-matrix, Eq. (4.8) with v J L , s , T , LST replaced 
by Tl'S't 1 lst-i an d the relation between the S- and T-matrices, Eq. (4.13), the amplitude 
induced by t> PC + t> PV can be expressed as 
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M s > m >t>,sm 8 t(E, 9) = V*^£ v / 2TTIe c y e a (L'(M S - M' S ),S'M' S | JM S ) 



JLL> 



(LO, SM S | JM s )Y L , (Ms _ M , )(6I) S *> aip) 5a '' a , (4.26) 
where again a(a') = LST(L'S'T'). 

C. Momentum- space formulation 

In order to consider the PC momentum-space Bonn (BONN) [11] and Nijmegen (NIJM- 
I) [10] potentials, it is useful to formulate the np scattering problem in p-space. One way to 
accomplish this is to solve for the if-matrix [33] 

K^Jp';p) =v J al Jp';p) 

+ 7fdttT4b';^4 a (^) > (4-27) 

where V denotes a principal- value integration, and the p-space matrix elements of the poten- 
tial are defined in Eq. (4.9). The integral equations (4.27) are discretized, and the resulting 
systems of linear equations are solved by direct numerical inversion. The principal-value 
integration is eliminated by a standard subtraction technique [34]. Once the K-matrices 
in the various channels have been determined, the corresponding (on-shell) ^-matrices are 
obtained from 



-i 



S J (p)= [l + 2i wK J (p;p)\ [l-2if,pK J (p;p)\ , (4.28) 

and from these the amplitudes M$'m' t',sm s t(E,0), Eq. (4.26), are constructed. 

Some of the studies of PV effects in the np system of interest here, specifically those 
relative to the np radiative capture, ^(7, n)p photo-disintegration at threshold and d(e, e')np 
electro-disintegration in quasi-elastic kinematics are more conveniently carried out in r- 
space, and therefore require r-space wave functions. To this end, one first re- writes Eq. (4.11) 
in a compact notation as 

W ^ ,P ^ = j(pr) - 2i fipj(pr) T J (p;p) 
r 

+ ^ rdkk 2 j(kr)^-T J (k;p) , (4.29) 
n Jo p 2 — k 

where the matrices [w J (r; p)] a ' iCC = w J a , a {r;p) and \j(pr)] a i jCe = Sa^ajuipr) have been intro- 
duced for ease of presentation. Then, by making use of the following relation between the 
off-shell T- and fT-matrices 

T J (p';p) = K\p'-p) -2ifi P K J (p';p)T J (p;p) , (4.30) 

which on-shell leads to 

T J (p;p) = \l + 2ifipK J (p;p)r 1 K J (p;p) , (4.31) 



13 



one can simply express the w J (r; p) matrix of solutions in terms of the previously determined 
fT-matrix as 



w J (r;p) 



x 



4/i f°° 2 ^7 

j(pr)-\ / dkk j(kr)— —K (k;p) 

n Jo p 2 — k 2 



1 + 2i up K J (p;p) 



(4.32) 



The Bessel transforms above are carried out numerically by Gaussian integration over a 
uniform p-grid extending up to momenta ~ 125 fm _1 . The computer programs have been 
successfully tested by comparing, for the PC Argonne v\$ (AV18) [9] and PV DDH [8] 
potentials, r-space wave functions as obtained from Eq. (4.32) and by direct solution of the 
Schrodinger equations, Eq. (4.15). 



D. The deuteron wave function 

The deuteron state has J=l and its normalized wave function is written in r-space as 
^d,m d (r) = E ^^STU LS T(r)yTs d ,j=i(r)Vo ■ (4-33) 

LST 

It has PC components with LST=010 and 210, the standard 3 Si and 3 D! waves (however, 
note again the unconventional phase factor i L , which makes the sign of the D-wave opposite 
to that of the S-wave), and PV components with LS'T=100 and 111, the x Pi and 3 P X waves 
(which are real functions because of the phase choice above). The radial functions are 
determined by solving the Schrodinger equation (4.15) in J=l channel with the boundary 
conditions Ulst(J') oc r L in the limit r — > and 

u ow (r) oc — , (4.34) 

,.. , (4-35) 

(4.36) 

in the asymptotic region. The asymptotic behavior of um(r) is identical to that of wioo( r ) 
above, and the constant k denotes the combination ^j2fi\Ed\, where \E^\ is the deuteron 
binding energy (2.225 MeV). 

In p-space the deuteron wave function is obtained from solutions of the homogeneous 
integral equations 

1 2 f°° 

U LS t(p) = -^ -TJTTy x - / dkk2 v LST,L'S'T'(P; k )u L 'S'T'(k) , (4.37) 

—P IK^) 7T JO L , S , T , 

and from these 



e -KT 




3 


oc 


1 + 




r 




KT 


-KT 


(" 


1 


^— 




r 




KT 



2 r°° 

ulst{t) = - dpp 2 j L (pr)u LST (p) ■ (4.38) 

7T JO 
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Figure 1 displays the functions ulst{^) obtained with the PC AV18 [9] (BONN [11]) and 
PV DDH [8] potentials (the values for the coupling constants and cutoff parameters in the 
DDH potential are those listed in Table I). The PC 3 Si and 3 Di components are not very 
sensitive to the input PC potential. For example, most of the difference between the AV18 
and BONN 3 Di waves is due to non-localities present in the one-pion-exchange (OPE) part 
of the BONN potential. In fact, it has been known for over two decades [35], and recently 
re-emphasized by Forest [36], that the local and non-local OPE interactions are related to 
each other by a unitary transformation. Therefore the differences between local and non- 
local OPE cannot be of any consequence for the prediction of observables, such as binding 
energies and electromagnetic form factors, provided, of course, that three-body interactions 
and/or two-body currents generated by the unitary transformation are also included — see 
Ref. [37] for a recent demonstration of this fact within the context of a calculation of the 
deuteron structure function A(q) and tensor observable T 20 (g) based on the local AV18 and 
non-local BONN potentials and associated (unitarily consistent) electromagnetic currents. 
This point was also stressed in Ref. [6]. 

The PV 3 Pi component is, in magnitude, much larger than the 1 Pi. This is easily 
understood, since the long-range 7r-exchange term in the DDH potential is non-vanishing 
only for transitions in which \T—T' | = 1, and therefore does not contribute in the x Pi channel. 
In this channel, however, the DDH p- and cu-exchange terms play a role. Note that, because 
of the short-range character of the associated dynamics, the AV18 and BONN 1 P 1 waves 
show considerably more model-dependence than the corresponding 3 Pi waves. 

Finally, in Fig. 1 the PV 3 Pi wave obtained with the AV18 and a truncated DDH 
potential, retaining only the short-range p- and cj-exchanges, is also shown. The comparison 
between the 3 Pi waves corresponding to the full and truncated DDH potentials demonstrates 
that this channel is indeed dominated by the 7r-exchange term in the DDH. 



V. PARITY- VIOLATING OBSERVABLES 

In this section we give explicit expressions for parity-violating (PV) observables in the 
np system, including the longitudinal asymmetry and spin rotation in np elastic scattering, 
the photon asymmetry in np radiative capture, and the asymmetries in deuteron photo- 
disintegration d(7, n)p in the threshold region and electro-disintegration d(e, e')np in quasi- 
elastic kinematics. 



A. Longitudinal asymmetry and spin rotation in np elastic scattering 

The differential cross section for scattering of a neutron with initial polarization m n is 
given by 

a mn (E,e) = E I « m ' P I M I m n™v) I 2 , (5-1) 

m P m' n m' p 

and the longitudinal asymmetry is defined as 

m6) - o + (E,9)+o„(E,9) ' (5 ' 2) 
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where ± denote the initial polarizations ±1/2. The total asymmetry A(E), integrated over 
the solid angle, then reads 



A(E) 



JdQa(E,9)A( y E, i 
JdQa{E,6) 



(5.3) 



where a = (a + + cr_) /2 is the spin- averaged differential cross section. The optical theorem 
allows A(E) to be simply expressed as 



A(E)=lmY: 

TT' L 



[(-: 



T+T' 



- 1 



Mut'ME, 0) + ^y— £ [l - (") 



T+T' 



\S+S' 



SS' 



X 



I m J2sm s t M S m s t,sm s t(E, 0) 
where in the equation above use has been made of the symmetry property 



Ms' t',sot(E, 0) 
(5.4) 



Ms',M' s ,T';S,Ms,t{E, 



:-) t+t \-) Ms - m ' s m s ,,_ m ,, t , s ,_ Ms , t {e^ 



(5.5) 



It is clear that the numerator of A(E) would vanish in the absence of PV interactions, since 
v PC , in contrast to t> PV , cannot change the total spin S or isospin T of the np pair. In 
particular, the long-range part of v PV due to pion exchange can only contribute to the first 
term in the numerator of Eq. (5.4), since it is diagonal in S, but non- vanishing for transitions 
\T-T | = 1. 

The transmission of a low-energy neutron beam through matter is described in terms of 
an index of refraction. A heuristic argument, outlined in Ref. [38], and the more rigorous — 
although less transparent — derivation presented in Ref. [33] show that a neutron with spin 
projection | m n ), after traversing a slab of width d of matter, is described by an asymptotic 
wave function given by 



e ip(z-d) e ipdn m „ 



(5.6) 



where p=p ra /2 is the initial relative momentum (assumed along the z- axis), and the index 
of refraction n mn is related to the density p of scattering centers in matter and the forward 
scattering amplitude. For the specific case under consideration here — neutron scattering 
from hydrogen — this relation reads: 



- 1 



2?rp 1 



V 2 2 



8=0 



Thus a neutron, initially polarized in the x-direction, 



(l+>- 



-»/V2 , 



(5.7) 



(5.8) 



having traversed a slab of matter, is described in the asymptotic region by a wave function 
given by 



e ip(z-d) e ipd(n + +n-)/2 



e ipd(n+-n-)/2 | _|_^ _|_ e ~rpd (n+-n_)/2 | j ^ . 



(5.9) 
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In the absence of v PV , the difference n + — n_ vanishes, since it is proportional to the sum 
over T,T' in the numerator of Eq. (5.4), while 



Im(n+ + p 



2p 



(5.10) 



where a is the spin-averaged cross section introduced above, and hence there will be an 
attenuation in the beam flux proportional to exp(— pdo). The real part of (n+ + n_)/2 
instead generates an un-observable phase factor. 

However, if PV interactions are present, then the real part of the (now non- vanishing) 
difference n + — n_ leads to a rotation of the neutron polarization by an angle explicitly 
given by [39-41] 



v 2 ^ 

f " m„ 



(+,m p | M | +,m p ) - (- 



M | -,m p ) 



(5.11) 



8=0 



B. Photon asymmetry in np radiative capture 

In the center-of-mass (CM) frame, the radiative transition amplitude between an initial 
continuum state with neutron and proton in spin-projection states m n and m p , respectively, 
and in relative momentum p, and a final deuteron state in spin-projection state rrid, recoiling 
with momentum — q, is given by 

JAm d , mnm >z> q) = (~q;"id I el(q) • j f (q) I pz,m n m p ) (+) , (5.12) 

where q is the momentum of the emitted photon and e^(q), A = ±1, are the spherical com- 
ponents of its polarization vector, and j(q) is the nuclear electromagnetic current operator. 
Note that p has been taken along the z-axis, the spin-quantization axis. 

The initial np continuum state, satisfying outgoing-wave boundary conditions, is related 
to that constructed in Sec. IV A via 

bz,m n m p )(+)=4=E(-) T+1 (^-^p \ SM s) \pz,SM s ,T)^ 

V ^ ST 1 1 

= v / 4^E^a V2LTl{ l -m n , l -m p \ SM s )(L0, SM S | JM s )\JM s ,a)^ , (5.13) 

Ja 11 

where in the first line the factor (— ) T+1 /v / 2 is from a Clebsch-Gordan coefficient combining 
the neutron and proton states to total isospin T, and in the second line e a = —(—) T e a and 
the states \JMj, a)^ have wave functions 

(r | J Mj , a) (+> = £ i" ^<^^, j( r) , (5.14) 

with a=LST and similarly for a'. The quantum numbers a and a' characterize the incoming 
and outgoing waves, respectively. 
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The CM differential cross-section for capture of a neutron with spin projection m n is 
then written as 



m p ,\,m d 

where 9 is the angle between z and q and 

a 



Jx. 



(+) 



m d ,m n m p 



(pz, q) 



(5.15) 



(TO = 



2n v 1 + q/m d 



(5.16) 



Here a is the fine-structure constant, is the deuteron mass, v is the relative velocity, 
v — p//i, and the photon energy q is given by 



q = m d 



-1 + 



i + 



\E d \ + 



m d 



2^ 



^ |£<zl + 



2/x 



(5.17) 



The photon asymmetry A 1 (9) is defined as in Eq. (5.2) with <r mn (9) replaced by cr^(9). 
By expanding the matrix elements of the current operator in terms of reduced matrix ele- 
ments (RMEs) of electric (E{) and magnetic (Mj) multipole operators as [28] 



(-q;m, | e* A (q) • jt(q) | JMj,a)W = -v^£(-i) 



'2Z + 1 



with 



£,(J,a) + AM,(J,a) 



X,(J,a) = (d, J=l || X, || J,a) (+) 



(JMj,^| lm^^-fl) 

(5.18) 



(5.19) 



and Xi—Ei or Mj, one finds, by retaining the 1 So and 3 Si channels in the sum over a in 
Eq. (5.13), the only relevant incoming waves in the energy regime of interest here (fractions 
of eV), 



A* (8) = a 'cos 9 , 



(5.20) 



where a 1 is given by 







i) + E* 1 ( 1 S ) 




+ Re 


E* 1 ( 3 S 1 )M 1 ( 3 S 1 ) 




M 1 ( 1 S )\* + 


E 1 ( 1 S )| 2 + |M 1 (3S 1 )| 2 + 


^l( 3 S!)| 2 



(5.21) 



and in Eq. (5.18) the d\ zV are standard rotation matrices [42]. The Ei(J,a) and Mi(J,a) 
RMEs should carry a superscript (+); it has been dropped for ease of presentantion. 

Several comments are now in order. Firstly, the photon asymmetry has the expected 
dependence on cos 9, since A 1 (9) oc <r n ■ q. Note that the contributions of higher order 
multipole operators with 1—2 have been ignored in the equation above. 

Secondly, because of the definition of the states in Eq. (5.14), a generic RME Xi(Ja) is 
expressed as 
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Xi(Ja) = ^X ; (Ja',a) , (5.22) 

a' 

namely as a sum over the contributions of outgoing channels ol corresponding to an incoming 
channel a; for example, 

MxOSq) = MxOSq, 1 S ) + Mi^Po, 1 S ) . (5.23) 

Thirdly, in the absence of parity-violating interactions, the only surviving RMEs are 
the M 1 ( 1 S , 1 S ), Mi( 3 Si, 3 Si), and M 1 ( 3 D 1 , 3 Si), and therefore the parameter a 1 vanishes. 
Furthermore, the Mi( 3 Si, 3 Si) RME also vanishes due to orthogonality of the initial and 
final states (in the limit in which isoscalar two-body currents are neglected), while the 
Mi( 3 Di, 3 Si) RME is suppressed in the energy regime of interest here. Thus, the standard 
result for the spin-averaged radiative capture cross-section, integrated over the solid angle, 
follows: 

^ = (4tt) 2 ( 7o|M 1 ( 1 So, 1 So)| 2 . (5.24) 

Lastly, when PV interactions are present, the analysis is more delicate, since then, in ad- 
dition to ad-mixing small opposite-parity components into the wave functions corresponding 
to v PC , these interactions also induce two-body terms in the electromagnetic current oper- 
ator, as discussed in Sec. III. Thus, j=j PC + j PV , where j PC includes the convection and 
spin-magnetization currents of single nucleons as well as the two-body currents associated 
with v PC , while j PV includes those terms generated by v PV . The multipole operators can then 
be written as Xu z —X^ + X Z PV , and those constructed from j PV have un-natural parities, 
namely (-)' for M[J and (- f +1 for Ef£ . Therefore, for example, the RMEs M PC ( 1 S , 1 S ) 
and M PV ( 1 So, 1 So) connect the PC x So state to, respectively, the PC and PV components 
of the deuteron. A straight-forward analysis then shows that, up to linear terms in effects 
induced by v PV either in the wave functions or currents, the parameter a 1 is given by 



-v^Re 


M^So/So) 


[E^Su 3 SO + £i( 3 Di, 3 Si) + E 1 ( 3 P 1 , 3 SO] 




M^So/So)! 2 



(5.25) 



where again terms containing the RMEs Mi( 3 Si, 3 SO and M 1 ( 3 D 1 , 3 Si) have been neglected. 
In the expression above, the RME £'i( 1 P 1 , 3 Si) has also been neglected, since transitions 
induced by the isoscalar electric dipole operator are strongly suppressed [43,44]. Thus, the 
only relevant transitions are those connecting the 3 Pi PV np state to the PC deuteron 
component and the 3 Si and 3 Di PC np states to the 3 Pi PV deuteron component. 

C. Helicity-dependent asymmetry in d(-j,n)p photo-disintegration 

The relevant matrix element in the photo-disintegration of a deuteron initially at rest in 
the laboratory is 

jii p ,Am d (P,q) = (_) (q;p,m„m p | e A (q) -j(q) | m d ) (5.26) 
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in the notation of Sec. VB above. Here | q; p, m n m p )^ represents an np scattering state 
with total momentum q and relative momentum p, satisfying incoming wave boundary con- 
ditions. Its (internal) wave function has the same partial wave expansion given in Eq. (4.10), 
except for the replacement w J a , a (r) — > [w^, a {r)\*. 

The cross section for absorption of a photon of helicity A, summed over the final states 
and averaged over the initial spin projections of the deuteron, reads: 

°l = E E ^ 6 (? + E d - E f ) ^ i j ( -! mp , Xmd (p, q) i 2 

"/'/'EE'.. I XMi-\j,a) + Et\j,a) \ 2 , (5.27) 



° Ja 1>1 



where Ef is the energy of the final state, 



q 2 p 2 



E f = — r + — , (5.28) 

f 2(m n + m p ) 2//' 1 ; 



and 



EE = -E E / r%- • ( 5 - 29 ) 

if ^ md mn m v J (271") 2 

Note that the factor 1/2 above is introduced to avoid double-counting the final states, and 
that the dependence upon the boundary condition of the continuum wave functions, the 
superscript (— ), has been reinserted in the RMEs of the electric and magnetic multipoles, 
namely 

X\'\j,a) = { -\j,a |] Xi || d, J = l) (5.30) 

and Xi=Ei or M h 

The resulting PV asymmetry, defined as P 7 = (cr+ — aV)/{a\ + al), is expressed as 



p7 _ 



J2ja J2l>l e a 


'M^~ ) {J,a)Ei~ ) *(J,a) + c.c.~ 






Mi~\j,a)\ 2 + 


Et\j,a) 


2 



(5.31) 



and therefore vanishes unless (i) the initial and/or final states do not have definite parity (as 
is the case here because of the presence of PV NN interactions) and/or (ii) the electric and 
magnetic multipole operators have unnatural parities (— ) l+1 and respectively, because 
of two-body PV electromagnetic currents associated with P V iVW interactions [7] . 

It is easily shown that in the inverse process p(n, j)d the expression for the photon circular 
polarization parameter is identical to that given above, but for the RMEs E^ and M/ - ** 
being replaced by the corresponding e\ + ^ and M.l + \ defined in the previous section. Indeed, 
by making use of the transformation properties of the states and electric and magnetic 
multipole operators under time reversal T, 

T\d,m d ) = (-) m ^ 1 \d,-m d ) , (5.32) 
T | J, Mj; = {-) M ^ J \ J, -Mj- , (5.33) 
TX Mz Tt = H 1 -*,,.,. , (5.34) 
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one finds the following relation for the RMEs: 

El +) (J,a) = (-) J+l El-\j,a), (5.35) 

and similarly for the Mj's. Hence, the circular polarizations measured in the direct and 
inverse processes are the same. 

D. Longitudinal asymmetry in d(e,e')np electro-disintegration 

The longitudinal asymmetry in the inclusive scattering of polarized electrons off a nuclear 
target results from the interference of amplitudes associated with photon and Z° exchanges 
as well as from the presence of parity-violating components in the nucleon-nucleon interac- 
tion. For completeness, we summarize below the relevant formulae. The initial and final 
electron (nucleus) four- momenta are labeled by /c M and k'^ {P^ and P' M ), respectively, while 
the four-momentum transfer is defined as = k^ — k ,fl = (u, q). The amplitudes for the 
7- and Z-exchange processes are then given by [45] 

M = -^(M 1 + M z ) , (5.36) 

M->=uYujl fi , (5.37) 

M Z = ^Z&u'-figP + , (5-38) 

where is the Fermi constant for muon decay, gy* = — 1 + 4sin 2 6 w and gff = 1 are the 
Standard Model values for the neutral-current couplings to the electron given in terms of 
the Weinberg angle 9w, u and u' are the initial and final electron spinors, and jjf and jff 
denote matrix elements of the electromagnetic and weak neutral currents, i.e. 

j}f=(f\P°(0)\i) = (p}MS fi (ci)), (5.39) 

and similarly for jjf . Here \i) and |/) represent the initial deuteron state and final np 
scattering state with incoming-wave boundary conditions (the (— ) solution), respectively. 
Note that in the amplitude M z the q 2 dependence of the Z° propagator has been ignored, 
since \q^\ <C m|. 

The parity-violating asymmetry A is given by the ratio of the difference over the sum 
of the inclusive cross sections dah/dVtdu for incident electrons with helicities h = ±1. It 
depends on the three-momentum and energy transfers, q and w, and scattering angle, 9 e , of 
the electron and conveniently expressed 

A = A 77 + A lZ • (5.40) 

Standard manipulations then lead to the following expression for the asymmetry in the 
extreme relativistic limit for the electron [31,45] 

v T 'R T , , . 

77 ~ v L RY f + v T R™ ' 1 ' 

, 1 G^g^v L Rr + gfv T R^ + g^v T ,R^ 
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where the i>'s are defined in terms of electron kinematical variables, 



v T = tan 2 (# e /2) + 
vr> = tan(0 e /2) 



2q 



2 ' 



\ 



tan 2 (fl e /2) + 



(5.43) 
(5.44) 

(5.45) 



The Rs are the nuclear electro-weak response functions, which depend on q and u, to be 
defined below. To this end, it is first convenient to separate the weak current j z ' a into its 
vector j°' a and axial-vector j 5,<T components, and to write correspondingly 



3fi = i°ff+i 5 ff = (P/i(q). J/i(q)) + (Pf<(q), j?i(q)) • 
The response functions can then be expressed as 

Rr(Q, ^ = £ £% + £ d - E f ) Re [p},(q)p^(q)" , 

* / 

Rr(q, W ) = EE% + ^- E f ) Re fc^q) • j^(q) 



i# b («Z, W ) = £ £ % + E d - Ej) Im [jj.(q) x jjj(q) 

* / 



(5.46) 

(5.47) 
(5.48) 
(5.49) 



where ^ is the ground-state energy of the deuteron (assumed at rest in the laboratory), 
Ef is the energy of the final scattering state, and in Eqs. (5.47) and (5.48) [Eq. (5.49)] the 
superscript a (b) is either 7 or (7 or 5). Note that there is a sum over the final states and 
an average over the initial spin projections of the deuteron. In the expressions above for 
the i?'s, it has been assumed that the three-momentum transfer q is along the z-axis, which 
defines the spin quantization axis for the nuclear states. 

The asymmetry induced by hadronic weak interactions, A 77 , is easily seen to be propor- 
tional to the interference of electric and magnetic multipole contributions as in Eq. (5.31), 
indeed 



(2tt) 2 « 



al-al 



(5.50) 



namely R^7 is related, of course for uj=q, to the difference of helicity-dependent photo- 
disintegration cross sections. Similar considerations to those in Sec. VC allow one to con- 
clude that this response would vanish in the absence of PV NN interactions (note, however, 
that in the present case there is, in addition to two-body PV currents, also a one-body PV 
term originating from radiative electro- weak corrections, the anapole current [23]). 



VI. CALCULATION 

In this section we briefly review the techniques used to calculate the PV observables in 
the np system — these are similar to those discussed most recently in Ref. [31]. 
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The deuteron wave function in Eq. (4.33) is written, for each spatial configuration r, as 
a vector in the spin-isospin space of the two nucleons, 

VWr) = X>il(r) |n), (6.1) 

n=l 

where \ n) = (p|)i (nT) 2 , (pfta • • • , (nJ-)i (pI)2 and are tne components of Vd,m d 
in this basis. The scattering wave function in Eq. (4.10) is first approximated by retaining 
PC and PV interaction effects in all channels up to a certain pre-selected J max and by using 
spherical Bessel function for channels with J > J max , and is then expanded, for any given 
r, in the same basis {| n)} defined above. The radial functions w a > tCe are obtained with the 
methods discussed in Sec. IV C. 

Matrix elements of the electromagnetic (and neutral- weak) current operators are written 
schematically as 

| O | ifc) = / dr^^^Om.nW^V) • (6-2) 

m,n 

The spin-isospin algebra is performed exactly with techniques similar to those developed in 
Ref. [17], while the r-space integrations are carried out efficiently by Gaussian quadratures. 
Note that no multipole expansion of the transition operators is required. 

Extensive and independent tests of the computer programs have been completed suc- 
cessfully. 



VII. RESULTS AND DISCUSSION 

In this section we present results for the longitudinal asymmetry and spin rotation in np 
elastic scattering, the photon asymmetry in the np radiative capture at thermal energies, and 
the asymmetries in the threshold photo-disintegration and quasi-elastic electro-disintegration 
of the deuteron. To provide an estimate for the model dependence of these results, we 
consider several different high-quality interactions fit to strong-interaction data, the Argonne 
v 18 (AV18) [9], Bonn (BONN) [11] and Nijmegen-I (NIJM-I) [10] interactions. We adopt 
the standard DDH [8] one-boson exchange model of the parity-violating (PV) interaction, 
and solve the Schrodinger equation for the scattering state and deuteron bound state with 
the methods discussed in Sec. IV. The values for the meson-nucleon coupling constants and 
cutoff parameters in the DDH model are those listed in Table I. Note that we have taken the 
linear combination of p- and u;-meson weak coupling constants corresponding to pp elastic 
scattering from an earlier analysis [6] of these experiments. The remaining couplings are the 
DDH "best value" estimates. As in the earlier analysis of pp scattering, the cutoff values in 
the meson-exchange interaction are taken from the BONN potential. 

It is also useful to introduce here some of the notation adopted in the following subsec- 
tions to denote variations on the DDH model defined above. The PV interaction denoted as 
DDHb corresponds to a DDH model with ir, p and uj weak coupling constants as specified 
by the "best value" set of Ref. [8], while the PV interaction denoted as DDH7T includes only 
the 7r-exchange term in the DDH model with the "best value" for the weak ttNN coupling 
constant. The remaining pion and vector- meson strong interaction coupling constants and 
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short-range cutoff parameters are as given in Table II for the DDHb model and are the same 
as in Table I for the DDH7T model. 

Finally, while the short-range contributions to the PV interaction should not be viewed 
as resulting solely from the exchange of single mesons, the six parameters of the DDH model 
are still useful in characterizing all the low-energy PV mixings. For example, two-pion 
exchange could play a role [46], however we assume that its effects can be included, at least 
at low energy, through the present combination of pion- and short-range terms. 

A. Longitudinal asymmetry and neutron spin rotation 

In this section we present results for the longitudinal asymmetry and neutron spin rota- 
tion in np elastic scattering. 

Figures 3-7 show the mixing parameters ef m induced by AV18 model in combination with 
the DDH, DDHb, and DDH7T interactions. Only those ef m induced by the PV interaction 
are displayed, namely for J=0 e\ 2 and for J > 1 ef m with 1—1,2 and m=3,4 in the notation 
of Tables III and IV. 

The definitions adopted for the phase-shifts and mixing parameters are those introduced 
in Sec. IV B. Up to linear terms in t> PV , the 5^ and e^ 1 values are not affected by weak 
interactions, and are determined solely by the strong interaction. They are identical to those 
listed in Ref. [9], but for two differences. Firstly, the Blatt-Biedenharn parameterization 
is used here for the S'-matrix [47] rather than the bar-phase parameterization of it [48] 
employed in Ref. [9] . Secondly, because of the phase choice in the potential components (see 
Eq. (4.16) and comment below it), the mixing parameters e^" 1 have opposite sign relative 
to those listed in Ref. [9]. 

The coupling between channels with the same pair isospin T is induced by the short- 
range part of the DDH interaction, associated with vector-meson exchanges; its long-range 
component, due to pion exchange, vanishes in this case. As a result, the mixing parameters 
in Figs. 3, 5, and 7, calculated with the DDH and DDHb models, are rather different, 
reflecting the large differences in the values for the some of the strong and weak coupling 
constants and short-range cutoffs between these two models, see Tables I and II. 

The mixing parameters between channels with |AT|=1, Figs. 4 and 6, in which the 
pion-exchange term is present, are still rather sensitive to the short-range behavior of the 
PV interaction, as reflected again by the differences in the DDH and DDHb predictions. 
However, this sensitivity is much reduced for the more peripheral waves, such as the 3 P 2 - 
3 D 2 and 3 F 2 - 3 D 2 channels. 

Figures 8-10 are meant to illustrate the sentitivity of the mixing angles to the input 
strong-interaction potential, which can be quite large, particularly in channels, such as the 
"Di^Pi. 

The total longitudinal asymmetry, defined in Eq. (5.3), is shown in Fig. 11 for a number of 
combinations of strong- and weak-interaction potentials. The asymmetries were calculated 
by retaining in the partial wave expansion for the amplitude, Eq. (4.26), all channels with 
J up to J max =6. There is very little sensitivity to the input strong-interaction potential. 
As also remarked in Ref. [6], this reduced sensitivity is undoutbly a consequence of the fact 
that present potentials are fitted to extended pp and pn databases with high accuracy. 
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Figure 12 shows that the total asymmetries obtained by including only the J=0 and 
1 channels (^o^Po, 3 Si- 3 Pi, 3 Di- 3 Pi, 'Si^Pi, and 3 D 1 - 1 P 1 ) and, in addition, the J=2 
channels, and finally all J channels up to J max =6. In the energy range (0-200) MeV the 
asymmetry is dominated by the J=0-2 contributions. 

For completeness, we present in Figs. 13 and 14 results for the angular distributions of 
the (PC) spin-averaged differential cross section and (PV) longitudinal asymmetry at center- 
of-mass energies of 20 MeV and 100 MeV. The asymmetry A(E,6) is defined in Eq. (5.2). 

The predictions for the neutron spin rotation per unit length, d0/d<i with defined in 
Eq. (5.11), are listed in Table V in the limit of vanishing incident neutron energy. The 
density of liquid hydrogen is taken as p=0.4 x 10 23 atoms-cm -3 . 

While results corresponding to different input strong interactions are within < 10% of 
each other, the calculated values show significant sensitivity to the short-range bahavior of 
the PV interaction, columns labeled DDH and DDHb. It is worth re-emphasizing that the 
longitudinal asymmetry in pp elastic scattering predicted by the DDHb model is at variance 
with that observed experimentally [6]. The short-range cutoff parameters and combinations 
of p- and cj-meson PV coupling constants in TT Z — 11, respectively h° p + h 1 + h 2 p / '\/6 and 
h-w + hi,, were constrained, in the DDH model, to reproduce this (measured) asymmetry [6]. 
An additional difference between the DDH and DDHb models is in the values adopted for the 
(PC) p- meson tensor coupling to the nucleon, 6.1 in the DDH (from the BONN interaction) 
and 3.7 in the DDHb (consistent with estimates from vector- meson dominance). Hence, 
the DDHb results are not realistic. Comparison between the DDH and DDH7T predictions, 
however, indicates that the neutron spin rotation is sensitive to the long-range part of v PV , 
and therefore a measurement of this observable would be useful in constraining the PV nNN 
coupling constant. 

Finally, there is a sign difference between the present results and those reported in 
Ref. [41]. It is not due to the different strong interaction potential used in that calcula- 
tion. Indeed, with the Paris potential [49] in combination with the DDHb model we obtain 
d0/dd=+8.88x 10~ 9 rad-cm _1 , the same magnitude but opposite sign than given in Ref. [41]. 

In order to understand this discrepancy, we have carried out a calculation of the neutron 
spin rotation, which ignores strong-interaction effects. It is equivalent to a first-order (in 
v PV ) perturbative estimate of this observable, and the corresponding results, listed in the last 
row of Table V (row labeled "plane waves"), demonstrate that strong-interaction distorsion 
effects are crucial, in fact they are responsible for flipping the sign of <fi. This is in contra- 
diction with the statement reported in the first paragraph after Eq. (6) of Ref. [41]: Avishai 
and Grange claim that the "plane-wave" prediction with the DDHb model is —6.0 x 10~ 9 
rad-cm -1 , namely it has the same sign as in their full calculation. 

The sign difference between the predictions obtained by either including or neglect- 
ing strong-interaction distorsion effects can easily be understood. For simplicity, con- 
sider the DDH7T model, in which case the relevant matrix element contributing to <fi is 
( 3 Px |w pv (DDHtt)| 3 Si), connecting the continuum T=0 3 Si and T=l Z Y> 1 channels. The es- 
sential difference between the un-distorted and distorted 3 Si wave functions is the presence 
of a node in the latter, thus ensuring its orthogonality to the deuteron 3 Si component. It is 
this node that causes the sign flip. 
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B. Photon asymmetry in np radiative capture at low energies 



The PV asymmetry a 7 in the 1 H(n,7) 2 H reaction at thermal neutron energies is cal- 
culated for the AV18, BONN and NIJM-I interactions. The asymmetry is expected to be 
constant for low-energy neutrons up to energies well beyond the 1-15 meV averaged in 
the experiment currently running at the LANSCE facility [50]. Each strong interaction 
model has associated two-body currents. For the AVI 8 we consider the currents from the 
momentum-independent terms — the 7r- and p-exchange currents from its Vq part — as well 
as from the momentum-dependent terms, as reviewed in Sec. III. Further discussion of the 
AV18 currents is given below. For the BONN and NIJM-I interactions, we retain only the 
7i- and p-exchange currents with cutoff parameters taken from the BONN model (A 7r =1.72 
GeV and A p =1.31 GeV), while we neglect contributions from other meson exchanges. In 
all calculations, however, the currents associated with the A excitation and uTvy transition 
have been included. 

The total cross section <r 7 is due to the well-known M x transition connecting the PC 
1 So np state to the PC deuteron state. The calculated values for each model are given 
in Table VI, both for one-body (impulse) currents alone and for the one- and two-body 
currents. In each case the largest two-body contribution, approximately two-thirds of the 
total, comes from the currents associated with pion exchange. The total cross section is in 
good agreement with experimental results, which are variously quoted as 334.2(0.5) mb [51] 
or 332.6(0.7) mb [52]. It would be possible to adjust, for example, the transition magnetic 
moment p 7 vA of the A-excitation current to precisely fit one of these values, here we simply 
choose a p 7 tva of 3 n.m., which is consistent with an analysis of j-N data at resonance. 

As discussed Sec. VB, the PV asymmetry arises from an interference between the M 1 
term above and the E x transition, connecting the 3 Pi PV np state to the PC deuteron 
state and the 3 Si PC np state to the 3 Pi PV deuteron state. The E\ transitions proceeding 
through the PV x Pi np or deuteron states are suppressed, because of an isospin selection rule 
forbidding isoscalar electric-dipole transitions and also because of spin-state orthogonality. 
In principle, there is a relativistic correction to the electric dipole operator, associated with 
the definition of the center-of-energy [44] . However, its contribution in transitions proceeding 
through the : Pi channel vanishes too, since the associated operator is diagonal in the pair 
spin. 

The calculated asymmetries are listed in Table VI. The results are consistent with 
earlier [53,54] and more recent [55,56] estimates, and are in agreement with each other at 
the few-per-cent level, which is also the magnitude of the contributions from the short-range 
terms. In particular, they show that this observable is very sensitive to the weak PV 7rNN 
coupling constant, while it is essentially unaffected by short-range contributions (in this 
context, see also Fig. 1). The E\ transition has been calculated in the long- wavelength 
approximation (LWA) with the Siegert form of the E\ operator (see, for example, Eq. (4.5) 
of Ref. [43]), thus eliminating many of the model dependencies and leaving only simple 
(long-range) matrix elements. In the notation of Sec. VB, the associated reduced matrix 
elements (RMEs) are explicitly given by 



£; 1 ( 3 S 1 , 3 S 1 ) + J E 1 ( 3 D 1 , 3 S 1 ) = i 



q 



2V67T 



f 

Jo 



w(r; 3 Pi 



w(r; 3 Si) 
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«;(r; 3 Dx) 
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(7.1) 

where the u>'s and it's denote the np continuum and deuteron radial wave functions defined, 
respectively, as in Sees. IV A and IV D (only the outgoing channel quantum numbers are 
displayed for the iu's). Corrections beyond the LWA terms in E\ transitions have been 
found to be quite small. For completeness, we also give the well known expression for the 
Mi RME, as calculated to leading order in q and in the limit in which only one-body currents 
are retained, 

M 1 ( 1 S , 1 S ) = i o i- (/i P -/i») rdrru^SxMrfSo) , (7.2) 
2V27rm Jo 

where the combination fi p — yU n =4.706 n.m. is the nucleon isovector magnetic moment. 

We have also calculated the E 1 contributions with the full current density operator j(x), 
namely by evaluating matrix elements of 

E 1X = - [ dxj(x) • V x 3l (qx)Y{l(±) , (7.3) 

where are standard vector spherical harmonics. To the extent that retardation correc- 
tions beyond the LWA of the E\ operator are negligible [43], this should produce identical 
results provided the current is exactly conserved. In order to satisfy current conservation, 
currents from both the strong (PC) and weak (PV) interactions are required, as discussed 
in Sec. III. In the following we keep only the 7r-exchange term in the DDH interaction 
(with their "best guess" for the weak tcNN coupling constant), and use the AV18 strong- 
interaction model. 

As reviewed in Sec. Ill, the PC two-body currents constructed from the v 6 part of the 
AV18 interaction (the ir- and p-exchange currents) exactly satisfy current conservation with 
it. The same holds true for the PV 7r-exchange currents derived from the DDH interac- 
tion in Sec. Ill A. However, the PC two-body currents originating from the isospin- and 
momentum-dependent terms of the AV18 are strictly not conserved (see below). The as- 
sociated contributions, while generally quite small, play here a crucial role because of the 
large cancellation between the (PC) v 6 currents from the AV18 and the (PV) 7r-currents 
from the DDH. This point is illustrated in Table VII. Note that the PC currents from 
A-excitation and ujh^i transition are transverse and therefore do not affect the E\ matrix 
element. However, they slightly reduce the PV asymmetry, since their contributions increase 
the Mi matrix element by ~ 1%. They are not listed in Table VII. 

The asymmetry is given by the sum of the two columns in Table VII, namely +0.17x 10~ 8 
(last row). This value should be compared to —5.02 x 1CT 8 , obtained with the Siegert form 
of the Ei operator for the same interactions (and currents for the Mi matrix element). As 
already mentioned, we have explicitly verified that retardation corrections in the E\ operator 
are too small to account for the difference. Thus the latter is to be ascribed to the lack of 
current conservation, originating from the isospin- and momentum-dependent terms of the 
AV18. 

To substantiate this claim, we have carried out a calculation based on a v$ reduction [57] 
of the AV18 (denoted as AV8), constrained to reproduce the binding energy of the deuteron 



- w(r; 3 Pi) 



u(r?Si) - -^u(r; 3 Di) 
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and the isoscalar combinations of the S- and P-wave phase shifts (note, however, that we do 
include in the AV8 the electromagnetic terms from the AV18, omitted in Ref. [57]). For the 
AV8 model, the n- and p-exchange currents from the isospin-dependent central, spin-spin, 
and tensor interaction components are constructed as for the AV18, and therefore are exactly 
conserved. However, the currents from the isospin-independent interaction, v b (rij)L ■ S, are 
derived by minimal substitution, 



Pi - ePi A(r; 



(7.4) 



where e and A are the electric charge and vector potential, respectively, and is the proton 
projection operator. The linear terms in A are written as — / dxj(x)-A(x), and the resulting 
spin-orbit current density — or, rather, its Fourier transform — reads 



j£g(q) 



v h ( rij ) 



(e iqr * Pi - e iq P 3 ) S 



x r. 



(7.5) 



In the case of the isospin-dependent terms, after symmetrizing v (r) [L • S , Tj • r,-] + /2, 
one obtains 



where 



^ (e iq ' r < Qj - v'^ Q, ) S x Tij , 



Tj + T z 



(7.6) 



(7.7) 



While minimal substitution ensures that the current is indeed conserved for the isospin- 
independent interaction, i.e. 



q-j£S(q)= v b ( rij )L ■ s , Pi (d) + Pj (d) 



(7.8) 



this prescription does not lead to a conserved current for the isospin-dependent one, since 
the commutator above generates an isovector term of the type 



i(T< X Tj) 



Jjz 



L S , e iq ri - e iq r 



(7.9) 



Physically, this corresponds to the fact that isospin-dependent interactions are associated 
with the exchange of charged particles, which an electromagnetic field can couple to. One 
can enforce current conservation by introducing an additional term [58], which in the case 
of jjj^ is taken as 



i(Ti X Tj); 



L • S , ri 



q • 



(7.10) 



+ 



The results obtained for the total cross section and PV asymmetry with the AV8 and 
pion-only DDH interactions and associated (exactly conserved) currents are listed in Ta- 
ble VIII. A few comments are in order. Firstly, the Mi cross section in impulse approxima- 
tion is ~ 30% smaller than predicted with the AV18 interaction. This is due to the fact that 
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the np singlet scattering length obtained with the AV8 (truncated) model is -19.74 fm, and 
so is about 15% smaller in magnitude than its physical value, -23.75 fm, reproduced by the 
AV18 within less than 0.1% [9]. 

Secondly, the enhancement of the Mi cross section in impulse approximation due to (PC) 
two-body currents, 9.3%, is essentially consistent with that predicted with the AV18. 

Lastly, the PV asymmetry obtained with the full currents is close to that calculated with 
the Siegert form of the Ei operator. The remaning ~ 1% difference is due to numerical 
inaccuracies as well as additional corrections from retardation terms and higher order mul- 
tipoles. Both of these effects are included in the full-current calculation. Note the crucial 
role played by the spin-orbit currents constructed above. 

C. Deuteron threshold disintegration with circularly polarized photons 

The photo-disintegration cross sections calculated with the AV18 and BONN models 
from threshold to 20 MeV photon energies are in excellent agreement with data [59]- [65], 
see Fig. 15. The model dependence between the AV18 and BONN results is negligible. In 
the calculations the final np states include interaction effects in all channels up to J=5 and 
spherical Bessel functions for J > 5, as discussed in the next section. 

In the energy regime of interest here, the (total) cross section is dominated by the 
contributions of E\ transitions connecting the deuteron to the np triplet P-waves. The 
Siegert form is used for the E x operator. Because of the way the calculations are carried out 
(see Sec. VI), it is conveniently implemented by making use of the following identity for the 
current density operator j(x), or rather its Fourier transform j(q), 



where in the first line the volume integral of j(x) has been re-expressed in terms of the 
divergence of the current, ignoring vanishing surface contributions, and in the second line 
use has been made of the continuity equation. Here p(x) is the charge density operator. In 
evaluating the matrix elements in Eq. (5.26) the commutator term reduces to 



where Pj is the proton projection operator introduced earlier, and relativisitc corrections to 
p(x), such as those associated with spin-orbit and pion-exchange contributions [43,44], have 
been neglected. 

We have also calculated the photo-disintegration cross section by using the expression 
given in Eq. (7.3) for the E\ operator, or equivalently by calculating matrix elements of the 
current j(q) without resorting to the identity in Eq. (7.11). The results obtained by including 
only the one-body terms and both the one- and two-body terms in j(q) are compared with 
those obtained in the Siegert-based calculation (as well as with data) in Fig. 16. The same 
conclusions as in the previous section remain valid here. Had the current been exactly 
conserved, then the Siegert-based and full j(q) calculations would have produced identical 




(7.11) 



dxx[# , p(x)] -> iq / dxxp(x) ~ ig J^Prj , 



(7.12) 
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results. The small differences in the case of the AV18 model, as an example, are to be 
ascribed to missing isovector currents associated with its momentum-dependent interaction 
components (see previous section). 

The PV photon polarization parameter P 7 , obtained with various combinations of PC 
and PV interactions, is displayed in Fig. 17, while its value at a photon energy ~ 1.3 keV 
above breakup threshold is listed in Table IX. All results presented below use the current 
operator in the form given on the right-hand-side of Eq. (7.11). Note that, as discussed in 
Sec. VC, the parameters P 7 for the direct d(j,n)p and inverse p{n^)d processes are the 
same. In the threshold region, a few keV above breakup, the expression for P 7 reduces to 

2Re[M 1 ( 1 So)P-( 1 S )] 
IM^Sq)! 2 

where, in the notation of the previous section, the M 1 ( 1 S ) RME is defined as in Eq. (5.23) 
and similarly for P 1 ( 1 So). In this energy region, the only relevant channel in the final np state 
has J=0, see discussion at the end of Sec. V B. Note that the combination of RMEs occurring 
in P 7 is different from that in a 7 , the photon angular asymmetry parameter measured 
in np radiative capture. Indeed, in contrast to a 7 , the photon polarization parameter is 
almost entirely determined by the short-range part of the DDH interaction, mediated by 
vector- meson exchanges (and having isoscalar and isotensor character [66]), see Table IX and 
Fig. 17. This is easily understood, since in the 1 So- 3 Po channel the pion-exchange component 
of the DDH interaction vanishes. Furthermore, the E\ transition connecting the 1 So np 
continuum state to the PV 3 Pi component of the deuteron, which is predominantly induced 
by the pion-exchange interaction, is strongly suppressed, to leading order, by spin-state 
othogonality. Higher order corrections, associated with retardation effects and relativistic 
contribution to the electric dipole operator, were estimated in Ref. [67] and were found to 
be of the order of a few % of the leading result arising from vector-meson exchanges. Some 
of these corrections are retained in the present study. 

The predictions in Table IX and in Fig. 17 display great sensitivity both to the strengths 
of the PV vector-meson couplings to the nucleon and to differences in the short-range 
behavior of the strong-interaction potentials, thus reinforcing the conclusion that these 
short-ranged meson couplings are not in themselves physical observables, rather the parity- 
violating mixings are the physically relevant parameters. 

Note that the ~ 5% decrease in P 7 values between the rows labelled "impulse" and "full" 
is due to the corresponding 5% enhancement of the M 1 transition connecting the PC 1 So 
and deuteron states, due to two-body terms in the electromagnetic current included in the 
"full" calculation. 

The results in Table IX are consistent both in sign and order of magnitude with those of 
ealier studies [68-71], remaining numerical differences are to be ascribed to different strong- 
and weak-interaction potentials adopted in these earlier works. Indeed, we have explicitely 
verified that by using the PC AV18 potential and the Cabibbo model for the PV potential [72] 
we obtain P 7 values close to those reported in Ref. [70,71]. However, our results seem to be 
at variance with those of Ref. [73] at photon energies a few MeV above above the breakup 
threshold. In particular, Table II in that paper suggests that at 10 and 20 MeV the dominant 
contribution to P 7 is from the PV pion-exchange interaction and that P 7 has the values 
—2.66 x 10~ 8 and —4.54 x 10~ 8 , respectively. This is in contrast to what reported in Fig. 17 
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(7.13) 



of the present work, curves labelled AV18+DDHb and AV18+DDH7T. There is a two-order 
of magnitude difference between the values referred to above and those obtained here. It 
is not obvious whether this is due to the use in Ref. [73] of the Hamada- Johnston [74] PC 
potential — a PV potential "close" to our model DDHb is adopted. 

The results in Table IX are consistent with the latest experimental determination, 
P 7 =(1.8 ± 1.8) x 1CT 7 [75], but about two orders of magnitude smaller than an earlier 
measurement [76]. 

Figure 18 shows the photon-polarization parameter obtained by including PV admixtures 
in the np continuum wave functions of all channels with J < J max and J max =0, 1, 2 and 5. 
In the energy range explored so far, P 1 is essentially given by the contributions of the J=0 
and 1 channels. 

Finally, Fig. 19 illustrates the effects of two-body terms in the electromagnetic current, 
written as in the right-hand-side of Eq. (7.11). The associated contributions are of the order 
of a few % relative to those from one-body terms. 

D. Deuteron electro-disintegration at quasi-elastic kinematics 

In this section we present results for the asymmetries A 77 and A lZ obtained by including 
one- and two-body terms in the electromagnetic and neutral-weak currents. Note, however, 
that only the PV two-body terms associated with 7r-exchange in the DDH interaction are 
considered in the present calculations (in addition, of course, to the PC terms discussed 
in Sec. III). The PV currents from p- and u;-exchange have been neglected, since they 
are expected to play a minor role due to their short-range character. One should also 
observe that at the higher momentum transfers of interest here, 100-300 MeV/c, relevant 
for the SAMPLE experiments [4,77], it is not possible to include the contributions of electric 
multipole operators through the Siegert theorem, these must be calculated explicitly from 
the full current. 

The A lZ contribution was recently studied in Ref. [31], where it was shown that two- 
body terms in the nuclear electromagnetic and weak neutral currents only produce (1-2)% 
corrections to the asymmetry due to the corresponding single- nucleon currents. The present 
study — a short account of which has been published in Ref. [7] — investigates the asymmetry 
originating from hadronic weak interactions. It updates and sharpens earlier predictions 
obtained in Refs. [78,79] — for example, these calculations did not include the effects of two- 
body currents induced by PV interactions. 

The present calculation proceeds as discussed in Sec. VI. We have used the AV18 or 
BONN models (and associated currents) in combination with the full DDH interaction (with 
coupling and cutoff values as given in Table I). The final state, labeled by the relative 
momentum p, pair spin and ^-projection SMs, and pair isospin T (Mr=0), is expanded in 
partial waves; PC and PV interaction effects are retained in all partial waves with J < 5, 
while spherical Bessel functions are employed for J > 5. In the quasi-elastic regime of 
interest here, it has been found that interaction effects are negligible for J > 5. 

In Figs. 20 and 21 we show, respectively, the inclusive cross section and the asymmetries 
A jZ and A J7 , obtained with the AV18 and DDH interactions, for one of the two SAMPLE 
kinematics, corresponding to a three-momentum transfer range between 176 MeV and 206 
MeV at the low and high ends of the spectrum in the scattered electron energy E', the 
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four-momentum transfer |g 2 | at the top of the quasi-elastic peak is ~ 0.039 GeV 2 . The rise 
in the cross section at the high end of the E' spectrum — the threshold region — is due to the 
Mi transition connecting the deuteron to the (quasi-bound) np 1 So state. Note that, because 
of the well-known destructive interference between the one-body current contributions orig- 
inating from the deuteron S- and D-wave components, two-body current contributions are 
relatively large in this threshold region. However, they only amount to a ~ 5% correction 
in the quasi-elastic peak region. 

In Fig. 21 the asymmetry A lZ labeled (l+2)-body — A lZ is defined in Eq. (5.42)— 
includes, in addition to one-body, two-body terms in the electromagnetic and neutral weak 
currents (in both the vector and axial- vector components of the latter). These two-body 
contributions are negligible over the whole E' spectrum. However, the (PC and PV) two- 
body electromagnetic currents play a relatively more significant role in the asymmetry Ay 7 , 
Eq. (5.41). 

In Fig. 22 we display separately, for the asymmetry A 77 , the contributions originating 
from i) the presence in the wave functions of opposite-parity components induced by the 
DDH interaction (solid curve) and ii) the anapole current and the PV two-body current 
associated with 7r-exchange (dashed curve). The latter are positive and fairly constant as 
function of E', while the former exhibit a pronounced dependence upon E' . Note that, up 
to linear terms in the effects induced by PV interactions, the asymmetry A 77 is obtained as 
the sum of these two contributions. 

The BONN model leads to predictions for the inclusive cross section and asymmetries, 
that are very close to those obtained with the AV18, as shown for A lZ and A 77 in Fig. 23. 
Thus the strong-interaction model dependence is negligible for these observables. 

In Fig. 24 we present results for the asymmetries corresponding to a four-momentum 
transfer \q 2 \ at the top of the quasi-elastic peak of about 0.094 GeV 2 , the three-momentum 
transfer values span the range (266-327) MeV over the E' spectrum shown. The calculations 
are based on the AV18 model and include one- and two-body currents. The asymmetry from 
7-Z interference scales with q 2 . and therefore is, in magnitude, about a factor of two larger 
than calculated in Fig. 21 where |g 2 | ~0.039 GeV 2 . The contributions to A 77 exhibit, as 
functions of E', a behavior qualitatively similar to that obtained at the lower |g 2 | value, see 
Fig. 22. 

In Fig. 25 we compare results for the contribution to A 77 due to the presence in the wave 
function of opposite-parity components induced by the full DDH and a truncated version 
of it, including only the pion-exchange term. Note that, up to linear terms in the effects 
produced by PV interactions, the other contribution to A 77 , namely that originating from 
(PV) one- and two-body currents, remains the same as in Fig. 22, since — as mentioned 
earlier — only the (PV) two-body currents associated with pion-exchange are considered in 
the present work. Figure 25 shows that the asymmetry A 77 is dominated by the long-range 
pion-exchange contribution. Hence, A 77 will scale essentially linearly with the PV irNN 
coupling constant. 

Finally, Fig. 26 is meant to illustrate the sensitivity of the asymmetry A 77 to those PC 
two-body currents derived from the momentum-dependent interaction components of the 
AV18 model (i.e., the spin-orbit, L 2 , and quadratic spin-orbit terms), see Sec. Ill and VII B. 
While these currents play a crucial role in the photon asymmetry in the np radiative capture 
at thermal neutron energy, they give negligible contributions to the present observable at 
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quasi-elastic kinematics. 

These results demonstrate that, in the kinematics of the SAMPLE experiments [4,77], 
the asymmetry from 7-Z interference is dominated by one-body currents, and that it is 
two-orders of magnitude larger than that associated with the PV hadronic weak interaction. 
Hence even the largest estimates of the weak irNN coupling constant will not affect extrac- 
tions of single-nucleon matrix elements. These conclusions corroborate those of the authors 
of Ref . [80] , who have carried out a similar study of the impact of hadronic weak interaction 
on quasi-elastic electro-deuteron scattering. 

VIII. CONCLUSIONS 

A systematic study of parity-violating (PV) observables in the np system, including the 
asymmetries in np radiative capture and ^(7, n)p photo-disintegration, the spin rotation and 
longitudinal asymmetry in np elastic scattering, and the asymmetry in electro-disintegration 
of the deuteron by polarized electrons at quasi-elastic kinematics, has been carried out by 
using a variety of latest-generation, strong-interaction potentials in combination with the 
DDH model of the PV hadronic weak interaction. We find that the model dependence of 
the np-capture asymmetry is quite small, at a level similar to the expected contributions of 
the short-range parts of the interaction. This process is in fact dominated by the long-range 
interaction components associated with pion exchange. A measurement of the np-capture 
asymmetry is then a clean probe of that physics. 

Similarly, we find that the asymmetry in the d(e, e')np reaction at quasi-elastic kinematics 
is a very clean probe of the electro-weak properties of individual nucleons. The processes 
associated with two nucleons, including PV admixtures in the deuteron and scattering wave 
functions and electromagnetic two-body currents induced by hadronic weak interactions, 
play a very small role at the values of momentum transfers explored so far [4,77]. 

We also find that the neutron spin rotation is sensitive to both the pion and vector-meson 
PV couplings to the nucleon, while exhibiting a modest model dependence, at the level of 
5-10%, due to the input strong-interaction potential adopted in the calculation. Thus a 
measurement of this observable [81], when combined with measurements of the asymmetries 
in np radiative capture [2] and pp elastic scattering [1], could provide useful constraints for 
some of these PV amplitudes. 

The asymmetry in the deuteron disintegration by circularly polarized photons from 
threshold up to 20 MeV energies is dominated by the short-range components of the DDH 
interaction. However, it also displays enhanced sensitivity to the short-range behavior in the 
strong-interaction potentials. Indeed, predictions for the asymmetry at threshold differ by 
almost a factor of two, depending on whether the Argonne v is or Bonn 2000 interactions are 
used in the calculations. Therefore, this observable cannot provide an unambiguous value 
of short-range weak meson nucleon couplings; however, they would be valuable in placing 
constraints on the hadronic weak mixing angles. 

Finally, the issue of electromagnetic current conservation in the presence of parity- 
conserving (PC) and PV potetials has been carefully investigated. In particular, in the 
case of the p(n,~f)d and d{^,n)p processes dramatic cancellations occur betweeen the con- 
tributions associated with the two-body currents induced, respectively, by the PC and PV 
potentials. 
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TABLES 



TABLE I. Values used for the strong- and weak-interaction coupling constants and short-range 
cutoff parameters of the 7T-, p-, and cj-meson to the nucleon in the DDH potential. 



all^ 




10 7 x h Q a 


10 7 x hi 10 7 x hi A a (GeV/c) 


TT 13.9 

p 0.84 
uj 20 




6.1 



-16.4 
3.23 


4.56 

-2.77 -13.7 
1.94 


1.72 
1.31 
1.50 


TABLE II. Values used for the strong- and weak-interaction coupling constants and short-range 
cutoff parameters of the 7T-, p-, and u;-meson to the nucleon in the DDH "best values" potential, 
from Ref. [8]. 


all^ 






10 7 x hi 


10 7 x h} a 10 7 x h 2 a 


A a (GeV/c) 


7T 13.9 

p 0.84 
u 20 




3.7 



-11.4 
-1.90 


4.56 

-0.19 -9.5 
-1.14 


2.4 
2.4 
2.4 








TABLE III. 


Labeling of channels. 




J 




1 




a 

2 3 


4 



1 
2 
3 




3 P 2 
3 D 3 




3p 

Po 

3 Di !Pi 
3 F 2 X D 2 
3 G 3 X F 3 


3 Pi 
3 D 2 

3 F 3 


TABLE IV. Classification of channel mixings for a given J: PC or PV if induced by w PC or 
i> PV , respectively. Note that no coupling is allowed between channels 3 and 4. 


J 


12 




13 


coupling 
14 


23 24 



1 

2 


PV 
PC 
PC 
PC 




PV 
PV 
PV 


PV 
PV 
PV 


PV PV 
PV PV 
PV PV 
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TABLE V. Neutron spin-rotation angle per unit length, in units of rad cm -1 , in the 

limit of vanishing incident neutron energy. Various combinations of strong- and weak interaction 
potentials are used. Also listed are the results obtained by ignoring strong-interaction effects, row 



labeled "plane waves" . 

DDH DDHvr DDHb 

AV18 IT09 K21 7A9 

NIJM-I 4.94 5.35 7.64 

BONN 4.63 5.18 7.35 

Plane waves -5.67 -6.87 -5.85 



TABLE VI. Total cross-section cr 7 and parity- violating asymmetry a 7 in the ftp radiative cap- 
ture at thermal neutron energies, obtained in various models. The asymmetry is reported for 
pion-exchange only (DDH7r) and full DDH (DDH) interactions. 



'(mb) 



x 7 x 10 8 



Interaction 
AV18 
NIJM-I 
BONN 



Impulse Current 
304.6 
305.4 
306.5 



Full Current 
334.2 
332.5 



DDHvr 

-4.98 
-5.11 
-4.97 



DDH 

-4.92 
-5.02 
-4.89 



TABLE VII. Cumulative contributions (in units of 10~ 8 ) to the PV asymmetry a 7 in the ftp 
radiative capture at thermal neutron energies for the AV18 interaction and pion-exchange-only 
DDH7T interaction. See text for explanation. 





AV18 (PC) Currents 


DDHvr (PV) Currents 


Impulse 


-15.3 




+VT 


-48.3 


44.2 


+P 


-40.4 


44.0 


+p-dependent 


-43.8 


44.0 
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TABLE VIII. Cumulative contributions to the total cross section <r 7 and PV asymmetry a 7 in 
the ftp radiative capture at thermal neutron energies for the AV8 and pion-exchange-only DDH7T 
interactions. Also listed is the asymmetry obtained with the Siegert form of the E\ operator. See 
text for explanation. 



a 7 (mb) a 1 x 10 ; 







AV8 (PC) Currents 


DDHvr (PV) Currents 


Total 


Impulse 


226.4 


-17.7 




-17.7 


+7T 


239.2 


-57.9 


51.3 


-6.60 


+P 


241.7 


-50.3 


51.1 


+0.790 


+SO 


241.7 


-57.0 


51.1 


-5.89 


+ A + W7T7 


247.4 


-56.3 


50.5 


-5.82 


Siegert E\ 








-5.76 



TABLE IX. Photon helicity-dependent asymmetries (in units of 10 -8 ) calculated with var- 
ious combinations of strong- and weak-interaction potentials at an incident photon energy of 
2.2259 MeV, about 1.3 keV above threshold. Predictions are listed obtained by including only 
one-body terms (impulse) and both one- and two-body terms (full) in the electromagnetic current, 
right-hand-side of Eq. (7.11). 





AV18+DDH (BONN+DDH) 


AV18+DDFLr 


AV18+DDHb 


Impulse 


5.44 (9.41) 


-0.035 


2.49 


Full 


5.19 (9.05) 


-0.037 


2.38 
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FIG. 1. The deuteron PC 3 Si and 3 Di and PV x Pi and 3 Pi radial wave functions obtained 
with the (PC) AV18 (black) or BONN (red) and (PV) DDH potentials. Also shown is the 3 Pi 
wave obtained with the AV18 and a truncated DDH potential (labeled "no tt in DDH"), including 
only p- and w-meson exchange contributions. For the phase convention, see text. 
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FIG. 2. Feynman diagram representation of the two-body currents associated with pion ex- 
change: solid lines, nucleons; dashed lines, pions; wavy lines, photons. Note that one interaction 
vertex is parity-conserving, while the other is parity- violating (PV). 
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FIG. 3. The ^o^Po mixing parameter obtained with the AV18 model in combination with 
either of two variations of the DDH model, labeled DDH and DDH "best values", see text. 
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FIG. 4. The 3 Si- 3 Pi and 3 Di- 3 Pi mixing parameters obtained with the AV18 model in com- 
bination with either of three variations of the DDH model, labeled DDH, DDH "best values", and 
DDH 7r-only, see text. 
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FIG. 5. The 3 Si- 1 Pi and 3 Di- 1 Pi mixing parameters obtained with the AV18 model in combi- 
nation with either of two variations of the DDH model, labeled DDH and DDH "best values" , see 
text. 
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FIG. 8. The 1 So- 3 Po mixing parameter obtained with the DDH model in combination with 
either the AV18 or BONN model. 
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FIG. 9. The 3 Si- 3 Pi and 3 Di- 3 Pi mixing parameters obtained with the DDH model in combi- 
nation with either of the AV18 or BONN model. 
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FIG. 11. The neutron asymmetry obtained with various combinations of strong- and weak 
interaction potentials, as function of the center-of-mass energy. 
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FIG. 12. Contributions to the neutron asymmetry obtained by including only the J=0 and 
J=l channels, and by adding the J=2 channels, and finally all J channels up to J m ax=6. The 
AV18+DDH potential combination is used, black solid line in Fig. 11. 
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FIG. 13. Angular distributions for the spin-averaged np (strong-interaction) cross section at 
center-of-mass energies of 20 MeV and 100 MeV, corresponding to the AV18 potential. 
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FIG. 14. Angular distributions for the neutron asymmetry at center-of-mass energies of 20 MeV 
and 100 MeV. The AV18+DDH potential combination is used, black solid line in Fig. 11. 
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FIG. 15. The deuteron photo-disintegration cross sections, calculated with the AV18 and 
BONN interactions, are compared to data. Note that the AV18 and BONN results are indis- 
tinguishable. 
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FIG. 16. The deuteron photo-disintegration cross sections, calculated with the AV18 and 
BONN interactions, are compared to data. Results obtained by including only one-body terms 
and both one- and two-body terms in the electromagnetic current are shown along with those 
calculated by using the right-hand-side of Eq. (7.11). The latter are the same as in Fig. 15. 




5 10 15 20 25 

co(MeV) 

FIG. 17. The photon helicity-dependent asymmetries obtained with various combinations 
of strong- and weak interaction potentials. Note that the predictions corresponding to the 
AV18+DDH7T potential combination are suppressed by roughly one order of magnitude relative 
to those corresponding to the AV18+DDH and AV18+DDHb models. All results are obtained by 
using the right-hand-side of Eq. (7.11). 
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FIG. 18. Contributions to the photon helicity-dependent asymmetry obtained by including PV 
admixtures in the wave functions of all channels up to J max , with J max =0, 1, 2, and 5. The 
AV18+DDH potential combination is used, black solid line in Fig. 17. Note that the curves labeled 
J = 2 and J = 5 are indistinguishable. 
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FIG. 19. The photon helicity-dependent asymmetries obtained with the AV18+DDH and 
BONN+DDH potential combinations by including only one-body terms and both one- and 
two-body terms in the right-hand-side of Eq. (7.11). 
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FIG. 20. The d(e, e')np inclusive cross section calculated, as function of the scattered electron 
energy E', with the AV18 interaction model. The electron incident energy is 117 MeV and its 
scattering angle e is 138.4°. Predictions are shown obtained with one-body terms alone and both 
one- and two-body terms in the electromagnetic current. 
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FIG. 21. The asymmetries and A 1 z calculated, as function of the scattered electron energy 
E' , with the (PC) AV18 and (PV) DDH interaction models. The other electron kinematical 
variables are as in Fig. 20. Predictions are shown obtained with one-body terms alone and both 
one- and two-body terms in the electromagnetic and neutral weak currents. 
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FIG. 22. Contributions to the asymmetry Ay 7 calculated, as function of the scattered electron 
energy E', with the (PC) AV18 and (PV) DDH interaction models. The other electron kinematical 
variables are as in Fig. 20. The solid line represents the results corresponding to the presence in the 
wave functions of opposite-parity components induced by the DDH interaction, while the dashed 
line represents the results due to the anapole current and the PV two-body current associated with 
7r-exchange. The total asymmetry shown by the solid line in Fig. 21 is obtained as the sum of 
these two contributions. 
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FIG. 23. The asymmetries Ay 7 and Ayz calculated, as function of the scattered electron energy 
E', with the DDH model in combination with either the AV18 or BONN model. The other electron 
kinematical variables are as in Fig. 20. Predictions are shown obtained by including one- and 
two-body terms in the electromagnetic and neutral weak currents. For the A^z asymmetry the 
AV18 and BONN calculated values are essentially indistinguishable. 
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FIG. 24. The asymmetry A^z and the two contributions to the asymmetry A 11 (notation as 
in Fig. 22) calculated, as function of the scattered electron energy E', with the (PC) AV18 and 
(PV) DDH interaction models. Predictions are shown obtained by including one- and two-body 
terms in the electromagnetic and neutral weak currents. Note that the electron incident energy is 
192 MeV and its scattering angle 9 e is 138.4°. 



64 



4 




60 70 80 90 100 110 

E'(MeV) 

FIG. 25. Contribution to the asymmetry A 77 , calculated as function of the scattered electron 
energy E', corresponding to the presence in the wave functions of opposite-parity components 
induced by either the full DDH or a truncated DDH model, consisting of its pion-exchange com- 
ponent only, in combination with the AV18 model. The other electron kinematical variables are as 
in Fig. 20. Predictions are shown obtained by including (PC and PV) one- and two-body terms in 
the electromagnetic current. 
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FIG. 26. The asymmetry Ay 7 calculated, as function of the scattered electron energy E', with 
the AV18 model in combination with a truncated DDH model consisting of its pion-exchange 
component only. The other electron kinematical variables are as in Fig. 20. Predictions are shown 
obtained by including the one-body terms alone and both the one- and two-body terms in the PC 
and PV components of the electromagnetic current (dashed and solid lines, respectively). Also 
shown are the results obtained by ignoring in the PC two-body currents those terms from the 
momentum-dependent components of the AV18 model (dashed-dotted line). 
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